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Introduction 


This  report  concentrates  on  progress  made  in  the  project  in  the  period  since  the  report  on  the 
first  phase  was  submitted  in  December  2004.  The  aims  of  the  project  were  extensively 
presented  in  the  first  report  and  are  only  repeated  briefly  here  since  the  goals  have  not 
changed.  The  aims  are:  (i)  To  use  data  acquired  in  the  DAWEX  campaign  of  2001  to 
investigate  in  detail  the  generation  of  gravity  waves  by  intense  convection  in  Northern 
Australia;  (ii)  To  use  equatorial  radar  observations  of  winds  in  the  MLT  (70-100  km)  to 
investigate  the  spectrum  of  gravity  waves  reaching  the  lower  thermosphere;  (iii)  To  make 
combined  observations  of  gravity  waves  and  ionospheric  irregularities  in  the  central  Pacific  to 
see  whether  gravity  wave  sources  can  be  identified  with  the  generation  of  specific  equatorial 
spread-F/plasma  bubble  events. 

In  the  first  report  the  focus  was  on  aims  (i)  and  (ii).  In  the  second  phase  the  focus  has  been  on 
completing  as  much  as  possible  the  work  on  aims  (i)  and  (ii)  and  making  a  transition  to  aim 
(iii),  which  is  the  focus  of  the  contract  awarded  under  this  project. 

(i)  Gravity  wave  studies  in  Northern  Australia 


A  series  of  campaigns  in  late  2001  used  radar,  radiosonde,  airglow  imager  and  modelling 
studies  of  gravity  wave  generation  by  strong  convection  in  the  vicinity  of  Darwin  (see  special 
issue  on  the  DAWEX  campaign  in  Journal  Of  Geophysical  Research,  Vol.  109,  D20,  2004). 
We  have  continued  analysis  and  interpretation  of  data  acquired  with  our  VHP  radar  on  the 
Tiwi  Islands,  and  in  particular  we  have  been  working  with  Dr  Joan  Alexander  (Colorado 
Research  Associates)  to  use  the  radar  observations  to  “calibrate”  her  high-resolution 
numerical  model  output.  This  has  proved  quite  successful  for  days  of  strong  convection 
(“Hectors”).  We  are  now  working  on  a  period  of  lower-level  convection  observed  during 
monsoon  conditions  (December  2001)  as  well  as  observations  made  with  support  from  this 
grant  during  the  dry  season  in  2004. 

We  are  also  comparing  model  predictions  of  wave  fluxes  in  the  lower  stratosphere  (20-35  km) 
with  observations  made  using  high-resolution  balloon  soundings.  An  unanticipated  problem 
has  been  the  impact  of  the  balloon  ascent  rate  (~5  ms'^)  on  the  retrieved  vertical  wavelengths 
of  waves  that  have  significant  vertical  phase  speeds.  We  are  assessing  the  effects  of  this 
“Doppler”  shifting  in  order  to  ensure  that  it  does  not  bias  our  estimates  of  wave  fluxes.  This 
work  is  being  carried  out  by  Peter  Eove  as  part  of  his  Ph  D  programme.  We  anticipate  one  or 
two  further  publications  to  eventuate  from  this  aspect  of  the  project. 

Finally,  it  is  worth  noting  that  outcomes  from  this  project  will  feed  into  another  major 
campaign  (TWP-ICE)  in  Northern  Australia  during  the  monsoon  season  in  January/February 
2006.  Together  with  Dr  Alexander  our  group  will  participate  with  the  aim  of  observing  and 
modelling  gravity  wave  generation  and  propagation  into  the  upper  atmosphere. 

(ii)  MF  radar  studies  of  gravity  waves  in  the  MET 

The  motivation  for  this  work  has  been  to  investigate  broader  aspects  of  the  gravity  wave 
source  spectrum  using  an  extensive  database  of  MF  radar  observations  at  several  near 
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equatorial  sites.  The  teehnique  is  based  on  a  Monte  Carlo  method  whereby  randomly 
generated  gravity  waves  with  known  speetral  properties  are  ray  traeed  through  the  middle 
atmosphere  and  the  resulting  synthesised  wave  varianees  eompared  with  observed  values.  The 
wave  speetrum  is  varied  until  satisfaetory  agreement  is  reaehed.  The  work  initially  deseribed 
in  the  phase  1  report  eoneentrated  on  two  equinoetial  months  (April  1993  and  April  1994)  as  a 
proof  of  eoneept.  The  teehnique  appears  to  be  robust  and  the  first  paper  has  been  aeeepted  for 
publieation  (Kovalam  et  ah,  2005).  One  reviewer  was  suffieiently  impressed  to  say  that  “I 
will  use  this  method  in  future.” 

During  the  seeond  phase  of  this  grant  we  have  extended  the  methodology  to  eover  other 
seasons.  The  results  are  similar  to  those  for  the  equinoetial  months  already  tested  and  show  an 
asymmetrie  speetrum,  as  illustrated  in  Figure  1.  The  main  feature  is  that  there  is  a  bias  toward 
waves  propagating  in  the  eastward  direetion.  While  the  exaet  reason  for  this  bias  is  not  known 
it  is  probably  due  to  a  shear  in  the  zonal  winds  in  the  lower  atmosphere.  We  are  now 
preparing  to  write  up  this  work,  ineluding  the  interaetions  between  gravity  waves  and 
atmospherie  tides  in  the  80-100  km  region.  While  the  latter  aspeet  is  ineidental  to  the  overall 
thrust  of  the  projeet  it  does  improve  our  understanding  of  the  role  of  gravity  waves  in  the 
MLT. 

(iii)  Extending  Gravity  Wave  Ray  Traeing  into  the  Thermosphere 

An  important  outeome  from  the  work  summarised  in  (ii)  is  that  waves  having  horizontal 
wavelengths  of  the  seale  apparently  assoeiated  with  equatorial  spread-F  do  appear  to  be 
present  in  all  seasons  in  the  lower  thermosphere.  The  next  stage  is  to  see  whether  these  waves 
ean  propagate  into  the  thermosphere  to  heights  of  at  least  200  km.  During  the  seeond  phase  of 
the  grant  we  have  been  extending  our  gravity  wave  ray  traeing  teehniques  to  inelude  the 
effeets  of  moleeular  viseosity  and  thermal  diffusivity.  The  approaeh  adopted  is  based  on  the 
work  of  Vadas  and  Fritts  (2004)  that  deseribes  a  new  method  of  ineluding  these  effeets  in  ray 
traeing. 

We  are  using  a  systematie  approaeh  in  developing  the  methodology.  Initially,  eonstant 
temperatures  were  used,  similar  to  the  examples  given  in  Vadas  and  Fritts  (2004).  The  next 
stage  has  been  to  use  more  realistie  temperature  profiles  derived  from  the  MSIS90  model.  The 
results  so  far  are  promising.  They  show  that  under  these  eonditions  only  longer  horizontal 
waves  with  large  vertieal  wavelengths  will  propagate  to  heights  of  200  km  or  above  (see 
Figure  2). 

The  next  stage  is  to  use  the  best  available  temperature  profiles  that  will  refleet  the  eonditions 
that  pertain  when  we  try  to  see  whether  speeifie  EsF/plasma  bubble  events  eould  be  triggered 
by  gravity  waves  propagating  from  known  eonveetive  sourees.  To  this  end  we  are  starting  to 
obtain  data  from  the  SABER  instrument  operating  on  the  TIMED  satellite.  It  will  also  be 
neeessary  to  use  the  best  available  wind  profiles,  sinee  the  baekground  winds  signifieantly 
affeet  the  speetrum  of  waves  reaehing  the  lower  thermosphere  and  subsequent  propagation  to 
higher  altitudes.  Several  different  ways  of  putting  together  suitable  wind  profiles  are  being 
investigated  in  the  phase  of  the  grant. 

Travel 

With  some  support  of  the  grant,  I  attended  the  C/NOFS  workshop  held  in  Estes  Park  in 
January  2005.  This  proved  invaluable  in  understanding  more  about  the  goals  of  the  mission 
and  the  role  we  might  play  in  supporting  it  onee  the  satellite  is  launehed.  Contaet  was  made 


with  Drs  Ed  Dewan  and  Keith  Groves  (AFRL,  Hanscom  AFB)  and  we  have  agreed  to 
collaborate  as  appropriate  in  support  of  C/NOFS. 


The  meeting  also  enabled  me  to  discuss  Christmas  Island  observations  of  EsF/bubbles  with 
Drs  Jonathon  Makela  (University  of  Illinois)  and  Roland  Tsunoda  (SRI).  They  have  airglow 
and  VHF  radar  instruments  directly  measuring  ionospheric  irregularities.  We  have  agreed  to 
collaborate  in  investigating  the  potential  for  gravity  wave  seeding  of  instabilities  by  pooling 
our  observations.  The  initial  focus  is  on  September  2004  when  spread-F  was  prevalent,  but 
not  ubiquitous,  at  Christmas  Island.  For  example,  strong  spread-F  was  observed  by  the  Cl 
VHF  radar  on  13  September  between  0530  and  1000  UT  (1900-0000  FT),  but  appeared  to  be 
absent  on  15  September.  As  part  of  the  3  phase  of  this  grant  we  have  started  to  investigate 
the  conditions  on  these  and  other  nights. 
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Figure  1:  Inferred  speetrum  of  convectively  generated  gravity  waves  in  the  20-120 
min  period  band  and  propagating  in  the  zonal  direetion  at  an  altitude  of  20  km  in  the 
vicinity  of  Christmas  Island  (2°N,  157  °W)  in  July  1993. 


Figure  2:  Height  profiles  of  four  gravity  waves  traced  through  the  middle  atmosphere 
into  the  thermosphere  in  a  no  wind  situation.  The  waves  have  horizontal  wavelengths 
of  100  km  (solid),  40  km  (dashed),  20  km  (dash-dot)  and  10  km  (dash-dot-dot-dot). 
The  panels  show  vertical  wavelength  (left),  normalized  wave  amplitudes  (center)  and 
normalized  body  force  (right).  The  key  result  is  that  only  the  100  km  wave  reaches 
heights  above  200  km  when  realistic  molecular  viscosity  and  diffusivity  is  included. 


